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Abstract. Fully differential cross sections are calculated for the ionization of H2 
by fast charged projectiles using a semiclassical model developed previously for the 
ionization of atoms. The method is tested in case of 4 keV electron and 6 MeV proton 
projectiles. The obtained results show good agreement with the available experimental 
data. Interference effects due to the two-center character of the target are also observed 
and analyzed. 
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In the last few decades there was a considerable development in the theoretical 



ON . description and experimental measurement of differential cross sections for charged 
particle impact ionization of atoms and molecules [IJ [2j [3j H] . Nowadays, the interest 



is focused on the in-detail analysis of the electron ejection from atomic or molecular 
targets [5j El U\ El El dD] • These analyzes may be performed by measuring and calculating 
fully differential cross sections which gives us the most complete information about an 
ionization process. These quantities describe the entire energy and angular distribution 
of the ionized electron, residual ion and projectile. 

Previously, based on the semiclassical impact parameter method we have 
constructed a theoretical model to calculate fully differential cross sections for single 
ionization of light atoms. This model takes into account the projectile-target nucleus 
interaction classically. The method was tested in case of the single ionization of helium 
produced by C 6+ ion projectile with an energy of E Q = 100 MeV/u and good agreement 
with the experiments was achieved in the scattering plane, while in the perpendicular 
plane a structure similar to that observed experimentally was obtained [TTl [12] . 

More complex and interesting features appear when the target is a molecule. Due 
to the multicenter character of the target, interference patterns occur in the ejected 
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electron spectra. This phenomenon was observed in the double differential cross section 
(DDCS) data for the ionization of hydrogen molecule by fast ions by Stolterfoht et al 
[T3] and has been analyzed theoretically by several groups pH [15j E3 E3, EEE] • 

Recently, fully differential cross sections were measured and calculated for the 
ionization of H 2 by both fast and slow electrons [9J [19], [20], and interference effects 
were analyzed in these data series. However, while several attempts have been made in 
the description of the interference effects in the double differential cross sections, and 
these effects are relatively well understood [MJ [16], it was shown recently [21], that an 
accurate description of the role of the projectile-target nucleus interaction remains a 
major challenge to theory. There are even much more details to clarify for these effects, 
if they appear in the total differential cross section. 

Interference effects due to indistinguishable diffraction of the incoming projectile 
from the two atomic centers were identified in the DDCS as a function of scattering 
angle for fixed ejected electron energy [22J. It was shown, that these kinds of interference 
structures can be more pronounced, that those in the ejected electron spectrum. 

The goal of the present paper is to adopt our previous semiclassical method [11] 
to calculate fully differential cross sections (FDCS) for the ionization of the hydrogen 
molecule, and to evidence the interference effects due to the two-center character of the 
target. In this model, the description of the molecule is similar to that from paper [23J, 
where total cross sections are calculated and are compared to the experimental data. 
Due to the limitations of a semiclassical model, we study the ionization induced by fast 
charged projectiles, and analyze the interference structures only in the ejected electron 
spectra. 

As a test case, the single ionization of the hydrogen molecule by 4 keV electron 
and 6 MeV proton projectiles is considered. FDCS are calculated and compared in the 
scattering plane with the experimental data of Cherid et al [21] and Dimopoulou et al 
[8]. Results in the perpendicular plane are also presented. In order to evidence the 
interference effects due to the two-center character of the target, the fully differential 
ionization cross sections for the H2 molecule and H atom are compared. 

2. Semiclassical theory 

In the semiclassical approximation the projectile is treated separately and it moves 
along a classical trajectory. This implies that only the electron system needs to be 
described by a time-dependent Schrodinger equation, while the projectile follows the 
classical laws of motion. 

2.1. General theory 

In order to study the ionization process of a small molecule induced by fast charged 
projectiles, first the ionization amplitudes have to be calculated. 

As described in [16], the first-order transition amplitude for a projectile with impact 
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parameter 6, velocity v and charge Z p , and a certain orientation of the molecular axis 
D may be written as 



a(6,k,D) = i-^ 

v 



+00 



dz e % 



^k(r) 



1 



|r-R 

where r is the position vector of the active electron and R = be 



>, (1) 
e z is the position 

of the projectile along the trajectory with z = vt. Here the origin is considered in 
the center of the molecule. The AE = Ei + is the energy transfer to the active 
electron, where Ei stands for the ionization energy. The initial and final state of the 
active electron is denoted by ^(r, D) and with k its momentum vector. In the 

calculations the Coulomb interaction 



-R 



is expanded into a multipole series. 



The transition probability for a given impact parameter and orientation of the 
molecular axis is the square of the transition amplitude modulus 

2 



w(b,k,~D) = a(6,k,D; 



(2) 



Because in experiments the molecular orientation is usually unknown, the transition 
probabilities are averaged over molecular orientation 



d \ dOo sin 60 w(b, k, D) 



(3) 



Finally, the triple differential cross section (TDCS) for the electron ejected in the 
energy range [E, E + dE] and into the solid angle dfl e and the projectile scattered into 
the solid angle dVt p can be expressed as 



d 3 a 



dVL e dVt p dE 



sin 6L 



db 



dQ r , 



w(b, k). 



(4) 



Because some experimental data [5] are differential relative to the perpendicular 
momentum transfer q^, we calculate also a FDCS expressed in terms of q_i_ and not 
of the scattering angle 



a* 



a 



dq±dk 



db 



dq 



w(b, k) . 



(5) 



An important part of the model is to assign impact parameter values to certain 
projectile scattering angles. As described in detail in [12], this may be achieved using 
the transverse momentum balance [2], which states that the momentum transfer is the 
sum of the transverse components of the electron's and residual ion's momenta. Further 
it is assumed, that the impact parameter is related to the momentum transfer to the 
residual ion, and we take into account the projectile-electron interaction separately. 
Thus, in this model, the projectile-target nucleus interaction is accounted for. By this 
point of view two main collision types are possible: (a) binary collision where most of the 
momentum transfer is taken by the electron and (b) recoil collision where most of the 
momentum transfer is taken by the target nucleus. Accordingly, two different impact 
parameter values have to be used in these collision regimes. The transition between 
these impact parameters is realized smoothly in the transition regions. 
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The projectile scattering is treated as a classical potential scattering problem in 
the field of the target system with nuclear charge Z t [25]. The simplest way to include 
the effect of the electrons around the target nucleus is to consider the potential to be a 
product of the Coulomb potential and the Bohr- type screening function [26] . Here it has 
to be mentioned, that using this potential we assume that the projectile scattering is 
produced by a spherical potential. This implies that in the description of the projectile 
scattering the two-center nature of the molecule is neglected. While in the analysis of 
the cross section as a function of the projectile scattering angle the correct description 
of the projectile-target nucleus interaction is important [21], we assume that the present 
approximation does not influence significantly the character of our results as a function 
of electron ejection angle. 



2.2. The particular case of hydrogen molecule target 

In the particular case of the hydrogen molecule the initial state is represented by a 
Heitler-London type molecular wavefunction. Written only for the active electron 

^(r, D) = N ( e -"l r - D / 2 l + e - Q l r+D / 2 l) , (6) 

where a is the effective charge, D is the vector associated to the internuclear distance 
and N denotes a normalization factor. In order to separate the dependence on the 
direction D of the molecular axis, this is expanded into a Legendre series [23]. In this 
series each term is characterized by the quantum numbers ^,mj. In the followings, the 
expansion coefficients depending only on the molecular distances will be denoted by 

The ^k( r ) continuum wavefunction is expanded into partial wave series, depending 
on angular momenta lf,rrif and phaseshift o~i f . The radial part Ri f (kr) is a wavefunction 
of the continuum electron moving in the mean field of the residual molecular ion 
approximated by the potential 

f I r > R 

V(r) = \ V r< l • (7) 

I D ' ' — 2 

This wavefunction is calculated numerically. 

By these considerations the ionization amplitude will be [23] 

a(6,k,D)=,— — x 

x ^ i~ lf e~ m f i f0 

1% v /(2/ / + l)(2/ c + l)(2/ J + l) lc0h ° 



E 3£U W^^W*)^ ^ b > D ) . (8) 



X 

Iflch 



where iVj and Nf are the normalization factors of the initial and final state target 
wavef unctions, respectively. The quantum numbers l c , m c corresponds to the multipole 
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expansion of the Coulomb interaction and 

G?l u (k,b,D)= f 00 dz e^ z Y lcmc (R) x 
j — 00 



-00 

/•OO ylc 

/ drr 2 R l} {kr)^ Cl .{r,D). (9) 
Jo r> 

is the integral over the projectile trajectory. 

Using this expression, the transition probability is calculated and averaged over 
molecular orientation. The FDCS for ionization of the hydrogen molecule will look as 

ofdcs = K — x 



N 



jHjr^ x 

X ifSefi ( 2l i + 1)^(2// + l)(2/' f + l)(2/ c +l)(2/^ + l) X 

r l f r l' f v ^/K+miJ^K+^J y nn N 

xF ;/(mc+mi) (fc)^ (m , +mi) (fc)G- b, D)G*™; h (k, 6, D) , 

where K is a constant depending on projectile charge, velocity and scattering angle, 
target charge and ejected electron momentum and it differs for the different kinds 
of FDCS mentioned earlier (see equations (j4j) and (jSJ)). While our results will be 
represented scaled to unity at their highest value in the scattering plane, in our following 
discussion the K constant will play the role of a simple scaling constant. 



3. Results and discussion 



In order to test the validity of the results of the semiclassical model for the H2 
ionization, calculated FDCS for fast electron and proton projectile impact are compared 
to the available experimental data [5| 12"^]. 

Fig. [1] shows fully differential cross sections for the ionization of the hydrogen 
molecule by 4.087 keV electron impact. The projectile scattering angle is 1° and the 
ejected electron energy is 20 eV. The top panel of the figure shows TDCS in the 
scattering plane, while the bottom panel shows the same data in the perpendicular 
plane. Experimental and theoretical data are each scaled to unity at their highest value 
in the scattering plane. Together with the results for the molecule (continuous line), 
TDCS are shown for the same collision process with atomic hydrogen target (dashed 
line). 

Further comment needs the fact that the TDCS in the scattering plane for atomic 
hydrogen target show shoulder structures at angles around 45° and 135°. Such structures 
cannot occur for pure first-order ionization from an isotropic Is state. As presented 
above, in order to take into account the projectile-target nucleus interaction, a certain 
impact parameter is assigned for every kinematic condition. Two different impact 
parameter values are used in the two collision regimes. The transition between these 



Semiclassical model for calculating ionization FDCS of H2 molecule 



6 



E = 4.087 keV;E =20eV;6 =1° 
p * • p 

= S 0° 90° 180° 270° 360° 




& »L ■ L ■ 4 ■ 4 ■ 1 

S, 90 180 270 360 

Electron ejection angle 9 e 



Figure 1. TDCS for ionization of H2 molecule in scattering (top) and perpendicular 
(bottom) planes. Experimental results of Chcrid ct al [24] are represented by points 
while our semiclassical results for the H2 molecule target arc drawn with continuous 
lines. TDCS values for the same collision process with atomic H target are represented 
by dashed lines. 

impact parameters is realized smoothly in the 6 e = 0°...50° and 8 e = 130°... 180° 
transition regions. The shoulders appear due to this impact parameter adaptation 
mechanism. 

In the scattering plane the theory reproduces well in shape the binary and the recoil 
peak structure. However, the semiclassical first-order results show a shift of the binary 
peak with some 10° relative to the experimental data, which have the angular resolution 
of less than 1° [51] . Such peak shifts indicates that the first order approximation may 
not be sufficiently accurate. 

Experimental data in the perpendicular plane for this particular process was not 
found. As in the case of our previous studies of the He ionization [12] . the semiclassical 
theoretical results show a double-lobe structure similar to that reported in the literature 
for slower collision processes [20] . 

In order to evidence the interference effects due to the two-center nature of the 
target, TDCS for the H2 molecule and H atom are compared in Fig. [TJ The only 
difference between these calculations is that the one-center wavefunction describing the 
H atom is replaced by a two-center wavefunction in case of the molecular target. At the 
first look the TDCS distributions for H 2 molecule and H atom are very similar. However, 
the recoil peak for the H2 target is larger than that for the H atom. The interference 
effects are analyzed later by the means of the interference factor. We note here, that 
we obtain the interference structures by comparing the theoretical results for H 2 and 
H. Because the agreement between our calculations and the experimental data is not 
perfect, we cannot say, that have identified the interference effects in the experimental 
data. 

We have identified the interference effects in the electron ejection spectrum, due 
to the coherent ejection from the two centers of the molecule. This effect is obtained 
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by describing the molecular electron by a two-center wavefunction. We assume, that 
neglecting the two-center character of the target in the description of the projectile 
scattering does not influence the obtained interference patterns. 

The interference structures in the ejected electron spectrum are analyzed through 
the interference factor defined by the ratio of the cross section obtained for the molecule 
and of two independent hydrogen atoms 0, [T7] 



°H 2 



2a H 

This interference effect may be detected in the TDCS distribution from non oriented 
molecules. The interference factor may be expressed by an analytical approximate 
formula [9] 

7 = 1 + 2^, (12) 

where q = q — k is the momentum imparted to the recoil ion. 

In Fig. [2] the interference factor obtained from our calculations (continuous line) is 
shown for the ionization process discussed earlier and it is compared to the analytical 
form ([12]) (dashed line). The top panel of the figure shows the interference factor in the 
scattering plane. As expected, it has oscillatory behavior with a strong maximum in the 
vicinity of the recoil peak. The semiclassical theory predicts another smaller maximum 
in the binary peak region, too, suggesting a more complicated behavior than that given 
by the analytical formula. In the perpendicular plane due to the kinematic conditions 
the q' momentum is constant: (1) the transverse component of the recoil momentum q' 
parallel to q is constant because the electron momentum is zero in that direction and 
q is fixed; (2) the transverse component of q' perpendicular to q is constant because 
q is zero in that direction and the electron energy is fixed. Accordingly, the analytical 
formula gives a constant interference factor. In contrast to this result, the semiclassical 
theory predicts some oscillations symmetric relative to the 180° direction (bottom panel 
of the Figure). 

TDCS distributions for other projectile energies and scattering angles are analyzed 
in Fig. [3J The top and middle panels present scattering plane TDCS results for 4.087 
keV electron projectile with 1.5° and 3° scattering angles. The ionized electron is ejected 
with 20 eV energy. The results are in good agreement with experiments both in the 
shape of the distribution and in the position of binary and recoil peaks. 

The bottom panel of Fig. [3] shows results for 4.167 keV projectile energy and a 
larger 8.2° scattering angle. The electron is ejected with a higher energy, of 100 eV. In 
this case the experiments show a narrow binary peak, which is well reproduced by our 
calculations. 

Fig. H] shows FDCS results of the semiclassical model for ionization of hydrogen 
molecule by 6 MeV proton projectile impact. In this case the experimental data 
of Dimopoulou et al [8] is differential in ejected electron momentum vector and the 
perpendicular momentum transfer vector. The comparison is made for different electron 
ejection energies and momentum transfers. 
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Figure 2. Interference factor for the process shown in Fig. Q] The continuous line 
represents results obtained by the present theory, while the dashed line shows the shape 
of the analytical form (fl2l) . 
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Figure 3. TDCS for ionization of H2 molecule in scattering plane for different 
projectile energies, scattering angles and electron ejection energies. The legend of 
this figure is the same as of Fig. [TJ 
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On one hand, for larger electron energies, of E e = 2.6 eV (top panels of the figure) 
the agreement between the experimental data and the present, semiclassical theory is 
reasonably good. Beside of this, it has to be noted that our first-order results show a 
shift of the binary peak with some 10°, but because of the experimental errors in the 
cross section this shift cannot be determined precisely. The cross section ratios for the 
binary and recoil peaks are also well reproduced. 

On the other hand, for lower electron energies, of E e = 0.2 eV (bottom panels 
of the figure) there is a discrepancy between the experimental and theoretical FDCS 
distributions. In experiments the cross section ratio of the binary and recoil peaks is 
almost 1. This feature is not reproduced by the semiclassical model. The explanation 
is the same as in the case of CDW-EIS model [8]: the FDCS distribution in the sub-eV 
region is influenced by the presence of the vibrational autoionizing channel (not included 
into the semiclassical model) which leads to the ejection of very low-energy electrons. 
The angular distribution of the autoionized electrons is essentially a dipolar one with 
respect to the momentum transfer axis. 

4. Conclusions 

In conclusion, the theoretical model based on the first order, semiclassical, impact 
parameter approximation used to calculate fully differential cross sections for single 
ionization of light atoms has been adopted to calculate fully differential cross sections 
for the ionization of diatomic molecules. The method has been tested in case of the 
ionization of H2 by 4 keV electron and 6 MeV proton projectiles. Except for some 
special cases, the obtained results in the scattering plane show good agreement with the 
experiments [H [24], and are in agreement with other theories [U El HSj. The double- 
lobe structure reported in the perpendicular plane is also reproduced. However, for low 
electron ejection energies the semiclassical model fails to reproduce the experimentally 
observed more symmetrical electron emission patterns in the scattering plane. In order 
to analyze the effect of the two scattering centers of the molecule on the cross section, we 
have compared the TDCS for the H2 with the TDCS obtained for independent atoms. 
Higher cross sections were obtained for molecules in the region of the recoil peak, and 
also the width of the binary peak was changed. These differences may be interpreted as 
interference effects. 

The aim of this work was to test that the semiclassical model, in its simple form, 
is able to treat complex colliding systems and describe interference patterns in the fully 
differential cross sections, too. 
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Figure 4. TDCS in scattering plane for the ionization of H2 molecule by 6MeV 
proton projectile for different momentum transfers and electron ejection energies. As 
in previous figures the continuous line represents the results of the semiclassical theory, 
while dots are representing the experimental data of Dimopoulou et al [5] . 
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